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We perform electronic structure and quantum transport studies of dangling bond loops created on H-passivated Si(lOO) surfaces 
and connected to carbon nanoribbon leads. We model loops with straight and zigzag topologies as well as with varying lenght 
with an efficient density-functional based tight-binding electronic structure approach (DFTB) . Varying the length of the loop or 
the lead coupling position we induce the drastic change in the transmission due to the electron interference. Depending if the 
constructive or destructive interference within the loop takes place we can noticeably change transport properties by few orders 
of magnitude. These results propose a way to engineer the closed electronically driven nanocircuits with high transport properties 
and exploit the interference effects in order to control them. 


1 Introduction 

During the last decades the miniaturization of electronic de¬ 
vices has steadily approached the atomic scale. As a result, 
a variety of new challenges arise, including the building of 
atomic scale circuit elements and logic elements on the ba¬ 
sis of individual molecules as well as gaining control over the 
inter-connects that ensure communication between the indi¬ 
vidual nanoscale components^HH. Hereby, scanning tunneling 
microscopy (STM) has revealed as a n invaluable tool to ma¬ 
nipulate matter at sub-nm lengthsEEl^ so that it has opened the 
possibility to engineer physical properties with atomic-scale 
precisionEl. In particular, a promising candidate for nanoelec- 
tronic applications are dangling bond wires (DBW), which are 
formed by selectively removing hydrogen atoms from Si(lOO) 
or Ge(lOO) passivated surfaces with the help of an STM tip; 
desorption of the H-atoms leaves behind an unpaired elec¬ 
tron (a dangling bond)El. This selective engineering at the 
atomic scale opens the fascinating possibility to design planar 
nanocircuits with complex geometry and tunable conduction 
propertiesCSEIl Additionally, DBW are inherently compati¬ 
ble with standard semiconductor based nanotechnologies, so 
that interfacing both domains may be a viable alternative in a 
near future. In real DB wires, however, the charge transport 
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efficiency may be strongly suppressed by Jahn-Teller distor¬ 
tions and buckling of the Si surface atoms^®^. Theoretical 
studies addressing the charge transport throu gh DB W mainly 
deal with charge transport along dimer rows^ ^h^^ i However, 
less attention has been paid to other, quasi 2D topol ogies, 
which have however been demonstrated experimentallyS^Uil 
Although electronic coupling perpendicular to the DB rows is 
smaller than along the dimer rows, it is non-zero and charge 
transport along DB loops may show interesting features like 
quantum interference effects. Motivated by this, we address in 
this investigation the problem of the charge transport proper¬ 
ties of realistic DB loops on Si(lOO) surfaces coupled to meso¬ 
scopic carbon nanoribbons acting as electrodes. We focus on 
intrinsic properties of the dangling-bond system like the loop 
topology and size and their influence on quantum interference 
effects, whose signatures show up in the quantum mechani¬ 
cal transmission function of the system. Since the DB loops 
need to be interfaced with the environment to probe charge 
transport, we also address the infiuence of the loop-electrode 
contact geometry and reveal a very sensitive dependence of 
charge transport on its local atomic structure. 

In our study of dangling-bond loops on H-passivated 
Si(lOO), we use a d ensity functional-based tight binding 
(DFTB) approachE^EUto perform structural relaxation as well 
as to compute the electronic structure of the different con¬ 
figurations dealt with. The DFTB approach allows to treat 
large structures (« 2000 atoms) with reasonable CPU time, 
while still keeping enough accuracy to provide a realistic in¬ 
sight into the properties of the system. From these calcula¬ 
tions we obtain the corresponding electronic density of states 
(DOS), the Fermi energy (Ep), and the Hamiltonian and Over¬ 
lap matrix elements. The latter two are then used as input 
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for a Green’s function based approach to compute the quan¬ 
tum mechanical transmission T{E) of the system. In brief, 
the standard Landauer-Biittiker formalisml22 is used, where 
the energy-dependent transmission function is computed as 
T{E) = rr(rL(^)G^(^)rR(^)G^(^)). Here, G^^^^ are re¬ 
tarded (advanced) Green’s functions of the system, rL,R = 

^ (^LR “ r) electrodes spectral densities, and ^l(r) 

are retarded (advanced) self-energy functions of the left L 
(right R) leads, which encode the electronic structure of the 
electrodes as well as the electronic coupling between electrode 
states and, in our case, the DB states. The self-energies are 
calculated with the well-known iterative Lopez-Sancho pro- 
cedureEI. 

We focus on two possible topologies of the DB loops, see[2 
for reference: (i) straight, when the Si surface is depassivated 
along the direction of the dimer row following a straight line, 
and (ii) zigzag, when the H atoms are removed in a zigzag 
way. The last one was predicted to have better charge trans¬ 
port properties in the case of infinite wiresl^l. There are two 
possible orientations of the loop on the surface, but in or¬ 
der to have a higher charge transport efficiency, the longer 
side of the loop needs to be oriented along the (100) direc¬ 
tion, since the electronic coupling perpendicular to the dimer 
rows is much smaller. During relaxation the Si DB surface 
atoms are buckled, showing a Peierls-like distortion, which 
may additionally reduce the transport efficiency when com¬ 
pared with the ideal unrelaxed loop. In a first step we ad¬ 
dress loops with a fixed linear dimension, consisting of five 
DBs along the dimer rows and four DBs in the perpendicular 
direction. The charge density associated to the DB states is 
localized very close to the Si surface as shown in the lower 
panel of^ To probe the transport properties of the DB loops, 
atomic-scale electrical contacts need to be engineered. One 
possibility is to use metallic gold pads^SES we have however 
chosen graphene nanoribbons as shown in due to their in¬ 
trinsic one-atom thickness and the variety of electronic fea¬ 
tures they display, which may add additional ways to tune 
electronic transport in hybrid carbon-silicon nanoscale de¬ 
vices. There exist already experimental studies with molec¬ 
ular nanowires and graphene nanoribbons on Si surfaces^SEll 
so that our proposed transport setup possesses potential exper¬ 
imental relevance. Due to their planar geometry and one-atom 
thickness, carbon nanoribbons can directly couple to single 
Si atoms on the DB loop, while remaining electronically de¬ 
coupled from the (passivated) Si substrate on which they are 
lyingESED. 

Apart from single atom contacts, it is also possible 
to have trapezoidal terminations of the graphene nanoribbon, 
i.e. more than one C atom can have a contact to Si atoms on 
the DB loop. This is also illustrated in panels b) and c) of|3. 
Naively, one may expect that increasing the number of atoms 
in the comer of the nanoribbon would improve the electronic 
coupling to the DB loop and thus increase the conductance (at 



Fig. 1 (a) Straight and (b) zigzag DB loops. The upper panels of (a) 
and (b) show the depassivated atoms highlighted in dark blue on an 
otherwise H-passivated Si(lOO) surface. The bottom panels of (a) 
and (b) display side views of the real part of electron wave function 
of the DB states in the Si band gap, where red color corresponds to 
the positive and blue - to the negative contributions. Although the 
largest spectral weight is on the dangling bond related surface states, 
there is still some contribution from deeper silicon layers. 


least at low energies). In fact, the local DOS is found to in¬ 
crease in the contact region (see SI), but at the same time the 
Fermi energy of the system is shifted and the spectral weight 
of the DB states in the Si band gap is modified, having a mix¬ 
ture of states coming from DBs and the leads. It turns then out 
that within a low-energy window ±0.2 eV around the Fermi 
level the leads ending with one C atom provide a more ef¬ 
ficient charge transport pathway than the two other cases with 
multiple contacts shown in|^ and[^. This behavior can be 
clearly seen in where we plot the transmission functions 
for the different contact geometries described above and con¬ 
sidered a straight loop consisting of five DBs along the dimer 
rows and four in the perpendicular direction. In the case (a) 
we obtain a high and narrow transmission peak near the Fermi 
energy close to the conduction band, while in the cases (b) 
and (c) quantum transport is strongly suppressed around the 
Fermi energy. The remaining peaks are related to the direct 
tunnelling between the DB rows and are discussed in SI. 

For the sake of clarity, we limit our following discussion to 
the case where a single carbon atom of the graphene nanorib¬ 
bon couples to a single atom on the Si DB loop, which cor¬ 
responds to [^. The nanoribbon slices for the Lopez-Sancho 
procedure were taken at a distance of few dimer rows away 
from the central part —including not only Si atoms, but also 
the first few atoms of the nanoribbon tip—, in order to mini¬ 
mize the influence on the DB states. This, together with peri- 
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Fig. 3 Transmission of a straight DB loop coupled to graphene 
nanoribbon electrodes terminated with (a) one (solid line in red, (b) 
two (dashed line in blue) and (c) three (dotted-dashed line in green) 
carbon atoms. 


odic boundary conditions, forces us to take quite a large sub¬ 
strate 8 Si atoms thick, 9 atoms width and 8 dimer rows long 
in the transport direction. 

Upon selectively removing the H atoms from the Si(lOO) 
surface to create the DB topologies, localized states emerge in 
the band gap of Si and the Fermi energy of the system (Ep) is 
shifted towards the band gap edge. It turns out that these local¬ 
ized states include contributions from both surface DB states 
and up to five layers of the Si substrate underneath Q; as a 
result, it is important to also include the substrate when com¬ 
puting the charge transport through the system^ZI. However, 
there is no direct tunneling between the loop arms oriented 
along the dimer rows, unless an electron reaches the connec¬ 
tion on top or bottom of the loop. The electronic DOS of the 
DB loop depends on the particular topology and it is different 
for straight and zigzag loops as shown in Eh and Eh- In the 
latter case more states, resulting from the symmetry breaking 
induced by the zig-zag conformation, appear on both sides of 
the Fermi energy within an energy window of roughly 0.2 eV. 
This eventually leads to a higher transmission (and current) 
for this topology, compare panels c and d in In the case 
of a straight loop the DB states are mainly created below the 
Fermi energy, in a smaller energy window (around 0.16 eV) 
closer to the conduction band. The corresponding transmis¬ 
sion displays therefore a very asymmetric behavior around the 
Fermi energy with strong suppression of the conductance al¬ 
ready few meV above the Fermi level (|^). 

Since DB loops form a close contour, one could expect to 
have a Mach-Zehnder like interferometer supported by the 
conducting surface states. In a semi-classical picture, elec¬ 
trons that enter the DB loop following different trajectories 
clockwise or counterclockwise have different phases when 
they reach the opposite lead. This may result in constructive 


or destructive interference effects, depending on the path the 
electron propagates. This effect might be expected only if the 
DB states take part in the transport, while tunneling between 
the loop side arms through the substrate will destroy any in¬ 
terference. Such quantum interference effects can be best re¬ 
vealed by varying the contact position between one of the elec¬ 
trodes and the DB loop. Onj^we show the effect of symmet¬ 
ric (solid line) and asymmetric (dashed line) coupling of the 
leads on the transmission function for straight (c) and zigzag 
(d) DB loops. In the case of asymmetric coupling, transport 
around the Fermi energy is dramatically suppressed for both 
topologies due to interference effects involving the DB elec¬ 
tronic states. The ratio between the transmissions for the sym¬ 
metric and asymmetric coupling geometries can be as large as 
10^. On the contrary, for energies farther away from the Fermi 
level the transmission is considerably less affected, since the 
corresponding electronic states have stronger bulk character. 

The interference effects can become more pronounced 
when tuning the electron phase by increasing the surface area 
of the DB loop. So far we were calculating electronic and 
transport properties of the full system that includes DB loop 
and part of the substrate as depicted on However, model¬ 
ing of larger loops becomes increasingly demanding from the 
computational point of view, so that instead of inverting large 
matrices to get the corresponding Green’s functions, we use 
numerically efficient recursive Green’s function techniquesEl. 
Details of the implementation are provided in the SI. 

Changing the length of the loop along the (100) direction is 
mapped then onto a change in the number of slices of Ny in the 
recursive Green’s function algorithm. Each slice is associated 
with one DB site and substrate beneath on each arm of the 
loop. Electrons moving clockwise or counterclockwise within 
the loop gain a different phase depending on the length. While 
increasing Ny we move the right lead along the dimer row ef¬ 
fectively ’’reading out” the electron phase difference along the 
loop. The influence of the loop length can be first illustrated 
with a simple toy model, see the top panel Of0 The model is 
described by a nearest-neighbor tight-binding Hamiltonian in¬ 
cluding only a single state per lattice site, but different hopping 
integrals between inequivalent sites: t\/t 2 = 2 .\j 2 /h = l .5^, 
and connected to square lattice semi-infinite leads. There are 
two clear effects in the transmission curve: even-odd oscilla¬ 
tions as a function of the loop length, and long-range oscil¬ 
lations that result from constructive (at A^ 3 ; = 16 and 91) and 
destructive (at = 51) interference processes. 

Going now back to realistic DB loops, we find for both 
straight and zigzag topologies the same qualitative behavior 
as discussed in the toy model, see the bottom panel Of0 The 
odd-even oscillations and long range interference effects are 
clearly seen if we fix the energy near a transmission resonance, 
as shown in panels a,b,d,and e in[^ while fixing it within a dip 
region (virtual tunnelling), we see an exponential decrease of 
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the transmission with increasing Ny, as illustrated in panels c 
and f of The transmission on panels a,b,d, and e may dif¬ 
fer by more than a factor of 10^ depending whether construc¬ 
tive or destructive interference takes place. Results concern¬ 
ing quantum interference signatures in the I-V characteristics 
of the DB loops with variable size can be found in the Supple¬ 
mentary Material. 

In summary, we have addressed from an atomistic per¬ 
spective charge transport signatures in dangling-bond loops 
on a Si(lOO) surface and contacted by mesoscopic graphene 
nanoribbons. Quantum interference effects have been demon¬ 
strated to be sensitively dependent on the connecting site of 
the electrodes to the DB loop. Loops built up with a zigzag ar¬ 
rangements of the depassivated atoms along the rows show a 
by far better low-energy (around the Fermi level) conductance 
when compared to the straight arrangement. Nevertheless, the 
strong asymmetry of the transmission around the Fermi level 
found in the latter case may have interesting consequences for 
the thermopower of the system, which is, in a first approxima¬ 
tion, determined on the behavior of the first derivative of the 
transmission around the Fermi level of the system. Addressing 
this issue requires however a separate study. Quantum inter¬ 
ference effects were also demonstrated to be only related to 
DB electronic states, so that their manipulation may allow for 
an additional control strategy of the electrical transport, since 
variations in the linear conductance of several orders of mag¬ 
nitude could be realized. This may be implemented by chang¬ 
ing the position of one of the leads, e.g. by using instead of a 
second graphene nanoribbon a functionalized STM tip, or by 
applying a voltage between different leads in a multi-terminal 
setup. Our investigation hints at the possibility of exploiting 
such effects in loop topologies for the design of planar-based 
atomic-scale electronics and its application in e.g. implement¬ 
ing logic gates. 

This work was partly funded by the EU within the projects 
Planar Atomic and Molecular Scale devices (PAMS, project 
nr. 610446). This work has also been partly supported by the 
German Research Foundation (DFG) within the Cluster of Ex¬ 
cellence ’’Center for Advancing Electronics Dresden”. Com¬ 
putational resources were provided by the ZIH at the Dresden 
University of Technology. 

Simulation details on different carbon nanoribbon leads 
coupled to the DB loops and technical details related to the re¬ 
cursive Green function approach used to model DBloops with 
variable size. 
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Fig. 2 Graphene nanoribbon electrodes terminated with one (a), two (b) and three (c) C atoms and corresponding electronic structure in an 
energy window of d=0.2 eV around the Fermi energy Ep . 



Fig. 4 Left panel: Electronic density of states (DOS) of (a) straight (red curve) and (b) zigzag (green curve) DB loops. Also shown as 
reference the corresponding DOS for the fully H-passivated Si(lOO) surface. Vertical dotted lines indicate the positions the Fermi level Ep in 
the different cases. Right panel: The corresponding quantum mechanical transmission functions at low energies, i.e. around Ep, for (c) 
straight and (d) zigzag DB loops. The dashed lines in both cases indicate the asymmetric connection of the DB loops to the graphene 
nanoribbons, see also the corresponding insets illustrating the different nanribbon-DB loop-nanoribbon contact configurations. 


6 I Journal Name, 2010, [vol]d-[To| 


This journal is ©The Royal Society of Chemistry [year] 

























































































































































N 

y 




Fig. 5 Top panel: Transmission function at energy Ep = 0.89 through the toy model loop with onsite energy equal to 4, measured in terms of 
hooping parameter between the sites of the semi-infinite squire lattice leads, connected to the loop, as a function of the loop length Ny. Bottom 
panel: Typical transmission functions through the straight (left) and zigzag (right) DB loops, as a function of the loop length Ny, calculated 
near the Fermi energy ?iiE — Ep = —0.02 (a), —0.027 (b) , —0.01 (c) — 0.05 (d), 0.05 (e), 0.0001 (f). 
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(a) (b) (c) 

Fig. 6 DOS the DB loop coupled to the graphene nanoribbon electrodes that terminate with one (in red, (a), two (in blue, (b)) and three (in 
green, (c)) C atoms. In black is the DOS of the fully passivated Si. The dotted lines show the corresponding position of the Fermi energy. 
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We provide here additional details on the graphene nanoribbons used as electrodes in the setups discussed in the manuscript. 
We focus on the local atomic structure of the electrode-DB contacts and their electronic DOS. We also shortly describe the 
recursive algorithm used to numerically solve the transport problem for larger loops. 

Multiple-contact topologies between graphene nanoribbons and DB atoms— The graphene nanoribbons were put on an H- 
passivated Si(lOO) surface at a distance of 3.5A. The nanoribbon edges are also passivated with H atoms, with exception of the 
triangular region, which mediates the electronic contact to the DB loop atoms. During the relaxation process, which ensures the 
maximal force component to be below 10“^ eV/A, the electrodes bend towards the closest Si dangling bond atom establishing 
a connection of the system to the leads with a 1.92A Si-C bond distance. This bond length is in good agreement with the 
experimental and simulation data reported before, see for example [Xu, Y.; He, K.; Schmucker, S.; Guo, Z.; Koepke, J.; Wood, 
J.; Lyding, J.; Aluru, N. Nano letters 2011, 11, 2735-2742]. 

The DOS for the different contact geometries described in the main text is shown in[^ We have considered a straight loop made 
of hve DBs along the dimer rows and four in the perpendicular direction. In the case of the single-atom contact, case (a), the 
position of the Fermi energy is only slightly shifted to +0.04 eV with respect to the DB loop alone, while in the cases (b) and 
(c) the relative change of the Fermi energy amounts to +0.26 eV. Moreover, for (b) and (c) the DOS around the Fermi level 
contains mainly contributions from states localized at the interface between the leads and the Si DBs (b,c). These states are 
weakly coupled to the rest of the DB loop, having an effective gap of 0.28 eV. These properties of the electronic structure are 
also reflected in the transmission function and are responsible for the strong suppression of quantum transport at low energies. 
These results reveal the way of how the lead should be connected to the DBs. The resulting states should be a mixture of both 
coming from the leads and DBs or having a small energy gap between them, such that they will be close to the Fermi level. In 
order to demonstrate role of direct electron tunnelling though the substrate we calculate the transmission function of the system 
made with two DBs on different dimer rows,coupled to the graphene nanoribbon leads. In this case transport of electrons can 
occur only by direct tunnelling perpendicularly to the dimer rows through the Si substrate. We compare it with the transmission 
function through the DB loop[^ One can see that around the Fermi energy transmission drops to which is much lower then 

in the both cases of destructive and constructive interferences. However around E — Ef = 0.21 the transmission peak, associated 
with the direct tunnelling between two DBs, is present in both cases. Exactly only this type of peaks remain in the case of charge 
localization on the interface between leads and DB loop discussed above. Recursive Green's function techniques— Modeling of 
larger DB loops is not practically achievable within a full hrst-principle framework, so that we use recursive techniques, building 
the system slice by slice in the transport direction as it is shown on[^ We construct the Green’s function separately for odd and 
even slices (between the dashed lines in[^ extracting the Hamiltonian and Overlap matrices from the hrst-principle calculations 
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Fig. 7 Transmission function around the Fermi energy for the transport within the DB loop (in red) and direct tunnelling between two DBs (in 
green). 


of the DB loop of smaller sizes. Here we used the assumption, supported by our DFTB calculations, that the Hamiltonian and 
Overlap matrices are approximately equal for each odd or even slice. This together with the computed self-energies and 
self- energies of bottom and top parts of the system Ztop,bottom allows us to calculate the transmission function between two leads 
fully taking into account the influence of the substrate as well as a varying number of slices Ny and loop geometry (straight or 
zigzag). 

We finally show inthe I-V characteristics of the DB loops with straight and zigzag topologies. The current is computed as 
I{V) = ^ /Too T(E) [/r(£’) — /l(£’)] dE, here T{E) is the zero bias transmission, and /l(R)(£’) are the left (right) electrode Fermi 
functions. As already predicted for the infinite wires, the current through the zigzag loop is considerably higher and sets in at 
lower bias (^ 0.05 V) when compared with the straight loop. The current mediated by the DB states dominates up to 0.2 V and 
0.15 V for zigzag and straight DB loops, respectively; afterwards charge transport through the conduction band starts and the 
current leaks into the Si substrate. 

The interference effects are less pronounced in the current vs. loop length plots shown in 10 since by integrating we average 
the phases over the bias energy window. However for small voltages, when the DB loops surface states start to conduct, for both 
loop topologies the quantum interference effects are still noticeable. 
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Fig. 8 Schematic illustration of the RGF procedure over Ny slices between two leads with self energies Each slice is defined as a part of 
the system between dashed lines, the solid line shows the part that contributes to the bottom- 



Fig. 9 Current as as function of bias voltage through the zigzag (green) and straight (red) DB loops, for a symmetric coupling of the graphene 
nanoribbons to the DB loop. 



Fig. 10 Current as a function of the DB loop length Ny, calculated at bias 0.12 V and 0.14 V for the straight (a) and bias 0.006 V and 0.06 for 
zigzag (b) topology. 
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